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(54) METHOD AND DEVICE FOR SUPPORTING COORDINATE GRID PREPARATION 

(57)Abstract: 

PURPOSE: To efficiently attain determining plural 
coordinate grids with different density by preparing a 
second coordinate grid whose interval is different from a 
first coordinate grid in the partial area of the first 
coordinate grid and preparing further new second 
coordinate grid with the display of corresponding 
relation. 

CONSTITUTION: Form data are inputted by using a 
general coordinate grid preparation supporting CAD 
system and the grid division of an area where analysis is 
executed in an analysis object. The end/not end of 
multiple grid preparation is judged and when it is judged 
to be not end, plural coordinate grid groups with 
different density are prepared and displayed. The second 
coordinate grid is prepared corresponding to the result 
of whether it is coarse or fine and it is stored in an 
external memory. Thus, the coordinate grid 
corresponding to the analysis object and an analysis 
purpose can be efficiently generated. 
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Partial English Translation of Japanese Laid-Open Patent 
Application No. H4-67 9 

It is assumed that a coordinate grid in a target analysis 
region is given by a general coordinate grid formation support 
CAD system or a known technique (for example, a curve coordinate 
grid is formed by coordinate transformation described in 
Japanese Patent Application No. 60-212985 and Japanese Patent 
Application No. 61-106466 ) . The subsequent procedures will be 
explained as a support method of the present invention. As a 
method of forming a multiple grid, a method of hierarchically 
forming a rough grid by thinning coordinate lines of a 
rough/dense grid was devised. As another method, the following 
method was devised. That is, the roughest grid is formed first, 
and coordinate lines are interpolated, so that a dense grid is 
hierarchically formed. In addition, a method of 

hierarchically the rough grid on the basis of a residual 
distribution of the rough grid and a method of forming a dense 
grid on the basis of a physical quantity distribution in the 
rough grid were devised. In particular, the later is 
effectively used when an optimum dense grid is formed while an 
initial value of the dense grid is calculated. In order to make 
it possible to form multiple grids, when a coordinate grid at 
a certain level is obtained in step 7, a menu "1. End, 2. Dense 
Rough, 3. Rough Dense, 4. Residual Display, 5. Physical 
Quantity Display" which is stored in an external memory 5b is 
displayed on the screen of a display device 2A by an arithmetic 
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unit 4a. An analyst selects a menu item from the menu to 
designate the menu item with a keyboard 6, so that the processes 
in steps 8 and 9 are executed, and the corresponding rough or 
dense grid is formed. 

A case in which an analyst forms a rough coordinate grid 
from a dense coordinate grid through many levels on the basis 
of her/his experience will be described below. A coordinate 
grid formed first is considered as a densest coordinate grid. 
In formation of a rough coordinate grid (step 8B) f the 
coordinate grid temporarily stored in an external memory 5a is 
displayed on the screen of the display device 2A again. It is 
assumed that this coordinate grid is called a coordinate grid 
(basic coordinate grid) at a current level. Rough grids are 
sequentially formed by thinning out several coordinate lines 
of the coordinate grid at the current level. FIG. 3 shows an 
example of a coordinate grid 14 at a current level displayed 
on the device of this embodiment. A procedure for forming a 
rough coordinate grid will be described below. An analyst , as 
shown in FIG. 4, designates one side (one coordinate line) 17 
of the coordinate grid 14 at the current level on the screen 
of the display device 2A by using the keyboard 6 or a light pen, 
designates a range (partial region) extending from 18 to 19 on 
the side, and picks the number of thinned coordinate lines in 
the range by a parameter input table 15 to designates the number 
of thinned coordinate lines. In the process in step 8B, the 
rough coordinate grid is calculated (by using, e.g., the 
technique described in Japanese Patent Application No. 61- 
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281773) in the partial region designated on the basis of the 
information obtained above. As shown in FIG. 4, a rough 
coordinate grid 16 obtained after the designated coordinate 
lines are thinned is displayed on the display devices 2A and 
2B. In this process , as shown in FIG. 4, a coordinate grid which 
includes the coordinate grid 14 at the current level and the 
rough coordinate grid 16 is obtained. On the other hand, as 
shown in FIG. 5, the thinned coordinate lines 14A (broken lines) 
and the rough coordinate grid 16 from which the coordinate lines 
14A are removed are displayed on the display device 2A. Another 
range on the same side or another coordinate line crossing the 
side and a range on the coordinate line are designated, and, 
similarly , a coordinate grid in which the designated coordinate 
lines are thinned is calculated and displayed. This operation 
is repeated until a coordinate grid which is satisfied on the 
basis of the experience of the analyst can be obtained, and a 
relation between the satisfied rough grid and a dense grid is 
stored in the external memory 5a, or the rough grid is also stored. 
The obtained rough grid is considered as a coordinate grid at 
a current level, and the above operation is repeated. In this 
manner, a combination (multiple grid) of rough and dense grids 
at multilevels can be formed. 

The designation of one side (or coordinate grid) is 
performed by picking one side with a light pen, a mouse, or the 
like. The designation of a range is performed by picking a 
coordinate line (coordinate surface in three-dimensional 
simulation) crossing the side 17. The number of thinned 
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coordinate lines is input by inputting the number of coordinate 
lines between thinned coordinate lines, by inputting the number 
of lines (surfaces) to be thinned, or by picking lines to be 
deleted . In the calculation of the rough coordinate grid , lines 
(surfaces) to be thinned are deleted, and grid point numbers 
are rearranged, and a relation between the rough coordinate grid 
and the dense coordinate grid is determined. In the display 
of the rough coordinate grid, the rough coordinate grid is 
emphatically displayed by colors, contour heavy lines , flicker, 
highlight, alternate display, parallel display, or the like. 
In the display of the relation, the rough and dense coordinate 
grids are displayed such that the coordinate grids are 
superposed (FIG.5). . 

When the analyst roughly divides a roughest coordinate 
grid on the basis of her/his experience to form a dense 
coordinate grid through many levels, a coordinate grid at a 
current level is the roughest coordinate grid (basic coordinate 
grid) at the present time. For example, it is assumed that the 
coordinate grid 14 in FIG. 6 is considered as a coordinate grid 
at a current level. In this case, according to a procedure for 
forming a dense coordinate grid, as shown in FIG. 6, one side 
17 of the coordinate grid 14 is designated, a range (partial 
region) extending from 18 to 19 on the side is designated, and 
the number of rough divisions in the range is designated. On 
the basis of the information, a dense coordinate grid (for 
example, by using the techniques described in Japanese Patent 
Application No. 61-106466, Japanese Unexamined Patent 
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Publication No. 62-263564, Japanese Patent Application 
No. 61-281773, Japanese Unexamined Patent Publication 
No. 63-136258, and Japanese Patent Application No. Hl-27462, 
and as shown in FIG. 6, the basic coordinate grid 16 and a dense 
coordinate grid 20 obtained by rough division performed by 
adding coordinate lines in a partial region are displayed on 
the display devices 2A and 2B. Another range on the same side 
17 or another coordinate grid crossing the side and a range on 
the coordinate line are designated, and, similarly, sub-divided 
coordinate grids are calculated and displayed on the display 
devices 2A and 2B. This operation is repeated until a 
coordinate grid which is satisfied on the basis of the 
experience of the analyst can be obtained, and the satisfied 
dense grid is considered as the densest coordinate grid. A 
relation of the rough coordinate grid to the dense coordinate 
grid and a fine coordinate grid are stored in the external memory 
5a. The coordinate data of the rough coordinate grid may be 
deleted. The densest coordinate grid is considered as a 
coordinate grid at a current level, and the above operation is 
repeated, so that the multiple grid can be formed. 

In the inputting of the number of subdivisions , the number 
of insertion divisions in a unit grid is input, or the number 
of all insertion divisions in a designated range is input. In 
the display of the rough coordinate grid, the rough coordinate 
grid is emphatically displayed by colors, contour heavy lines, 
flicker, highlight, alternate display , parallel display, or the 
like on the display devices 2A and 2B. In the display of the 
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relation, the rough and dense coordinate grids are displayed 
on the display devices 2A and 2B such that the coordinate grids 
are superposed (FIG, 6). 

On the basis of physical quantity calculation using the 
coordinate grid at the current level, a residual distribution 
for the coordinate grid is obtained (step 9B). When a rough 
coordinate grid is formed on the basis of the residual 
distribution (step 8B), a rough coordinate grid is set by a 
method (will be described below) such that a residual 
distribution in a control volume for one unit grid of the rough 
coordinate grid has a predetermined positive value (residuals 
have positive values) or a predetermined negative value 
(residuals have negative values) as much as possible. For 
example, when the numbers of vertical and horizontal unit grids 
included in a negative region (region in which the residuals 
have positive values) 21 and a negative region (region in which 
the residuals have negative values) 22 are equal to each other 
(FIG. 7A) , the unit grids are integrated in units of the positive 
region 21 and the negative region 22 adjacent to each other to 
a rough coordinate grid as shown in FIG. 7B. In the coordinate 
grid at the current level, when the numbers of unit grids 
included in the negative region 21 and the negative region 22 
of the residuals are considerably different from each other in 
the coordinate grid at the current level (FIG. 7C), the unit 
grids of the coordinate grid at the current level are integrated 
in units of the regions to form a rough coordinate grid as shown 
in FIG. 7D. More specifically, when the unit grids of the 
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coordinate grid at the current level are integrated with each 
other in units of the positive region (or negative region) of 
the residuals, one unit grid in the rough coordinate grid is 
formed. The coordinate lines of the rough coordinate grid are 
located on a boundary between the negative region 21 and the 
negative region 22 as much as possible. 

The sizes of the residuals are not classified into the 
positive region and the negative region, and the sizes of the 
residuals can also be classified into three or more regions in 
units of predetermined level ranges. In this case, the unit 
grids of a coordinate grid at a current level are integrated 
with each other in units of level ranges, and a rough coordinate 
grid is formed by the same manner as described above. 

A procedure for forming a multiple grid in consideration 
of a residual distribution is as follows. A case in which a 
residual distribution is classified into the positive region 
21 and the negative region 2 2 will be exemplified. A physical 
quantity on a coordinate grid at a current level is set, the 
physical quantity is repeatedly calculated, and the repetition 
is ended by an appropriate number of times (step ?B). The 
obtained residual distribution (the positive region 21 and the 
negative region 22 ) and the coordinate grid at the current level 
which is used in the calculation as shown in FIG. 8 are displayed 
on the display devices 2A and 2B (step 9C). 

In step 8B, retrieval for checking whether a rough 
coordinate grid corresponding to the displayed residual 
distribution is stored in the external memory 5a or not is 
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performed. When the corresponding rough coordinate grid is 
stored, the rough coordinate grid is called , and the rough 
coordinate grid is emphatically displayed on the superposing 
display of the residual distribution and the coordinate grid 
at the current level . When the corresponding rough coordinate 
grid is not stored, in step 8B, it is displayed on the display 
device 2A that the rough coordinate grid is not stored. An 
analysts designates coordinate lines to be deleted with 
reference to the display data of the residual distribution and 
the coordinate grid at the current level by the same manner as 
described above. In response to the designation information, 
in step 8B, a rough coordinate grid is formed from the dense 
coordinate grid by the thinning procedure, and, similarly, the 
rough coordinate grid is emphatically displayed (FIG. 9). 
Coordinate lines 14B indicated by chain lines in FIG. 9 are 
thinned lines. As described above, according to this 
embodiment, since the sizes of the residuals can be reflected 
on formation of a rough coordinate grid, a rough coordinate grid 
(multiple grid) having unit grids having different sizes can 
be easily formed as shown in FIG. 9. Since a residual 
distribution is displayed, coordinate lines to be thinned in 
the formation of the multiple grid can be easily known by an 
analyst . The analyst decides on the basis of her/his experience 
whether the emphatically displayed rough coordinate grid is 
appropriate or not (FIG. 9 shows an appropriate rough coordinate 
grid). When the obtained rough coordinate grid is not 
appropriate, the analyst designates coordinate lines to be 
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thinned from the coordinate grid at the current level again, 
and a rough coordinate grid is formed by the thinning procedure 
in step 8B. The rough coordinate grid is corrected by an 
operation (e.g., using the technique described in Japanese 
Patent Application No. 61-281773 and Japanese Patent 
Application No. Hl-27462) such as an operation of deleting or 
adding a coordinate line on the obtained rough coordinate grid. 
The formed rough coordinate grid is emphatically displayed on 
the display devices 2A and 2B such that the rough coordinate 
grid is superposed on the residual distribution and the 
coordinate grid at the current level (FIG. 9) . When the analyst 
inputs a signal representing that the rough coordinate grid is 
appropriate, a relation of the rough coordinate grid to the 
dense coordinate grid and the rough coordinate grid are stored 
in the external memory 5a. FIG. 10 shows a residual 
distribution obtained when a Poisson equation is solved on a 
13x13 coordinate grids 14C at a current level and the coordinate 
grids 14C at the current level used in the calculation. The 
residual distribution is indicated by contour lines 30. FIG. 
11 shows a rough coordinate grid (solid line) 16A formed by the 
above method together with the residual distribution in FIG. 
10 and the coordinate grid (chain line) 14C at the current level. 

As described above, upon completion of the formation of 
the multiple grid in consideration of the residual distribution, 
when YES is decided in step 7C, the process in step 10 is executed 
by using the obtained multiple grid. In the display example 
in FIG. 9, the rough coordinate grid 16 having an uneven density 
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includes 11 unit grids having different sizes depending on the 
residual distribution. Simultaneous linear equations A 1 x 1 = 
b* are formed in each of the unit grids, and the equations are 
repeatedly solved until the residuals have predetermined values 
or less. As an initial value of the equations, a value obtained 
by solving predetermined simultaneous linear equations Ax = b 
to the unit grids of the coordinate grid at the current level 
is used. In this embodiment, the rough coordinate grid having 
unit grids having different sizes can be formed depending on 
a residual distribution. For this reason, when the 
simultaneous linear equation A' x 1 = b 1 are solved, the influence 
of the residuals can be eliminated within a short time. 
Therefore, convergence time is short, and an obtained physical 
quantity can be calculated at high accuracy. 

In the display of the residual distribution, the 
calculation grid, and the rough grid, a superposing display and 
a parallel display can be employed. 

When a further dense coordinate grid is formed on the 
basis of the calculated physical quantity, the physical 
quantity used in the calculation and the coordinate grid at the 
current level are displayed such that the physical quantity is 
superposed on the coordinate grid. A further dense coordinate 
grid is formed from the coordinate grid at the current level 
used in the calculation is formed by the subdividing procedure 
in step 8D in FIG. 1. At this time, when coordinate line 
positions of the coordinate grid used in the calculation are 
not equal to coordinate line positions of the dense coordinate 
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grid, and when employment of the dense coordinate grid is 
designated by the analyst, by using the dense coordinate grid 
as the coordinate grid at the current level, the coordinate grid 
used in the calculation and all coordinate grids which are 
rougher than the coordinate grid used in the calculation are 
formed again (process in step 8D). 
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